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1.0 


Assumptions 


In  this  study  I generally  made  the  following  assumptions  and 
restrictions  on  what  was  considered: 


We  would  have  available  only  a single  three- 
component  seismometer; 

We  wish  to  find  the  frequency  domain  amplitude 
of  both  Love  and  Rayleigh  waves  of  an  explosive 
signal; 

The  origin  time  and  great  circle  azimuth  to 
the  explosive  source  is  known; 

Interference  may  be  natural  or  an  interfering 
eartl  quake; 

The  period  band  of  20  to  60  seconds  is  of 
interest,  with  the  primary  interest  in  40-second 
periods. 


Literati!  re 


2.  0 


2.  1 The  Nature  of  Interference 

A.  The  natural  background  (excludes  identifiable 
interfering  earthquakes): 

1.  At  the  shorter  period  end  (20  to  30  seconds), 
the  natural  noise  has  two  major  components.  Approximately  half 
the  noise  power  comes  from  propagating  fundamental  Rayleigh 
waves  generated  by  ocean  wave  action.  The  level  of  this 
noise  is  variable  with  meteorological  conditions  in  the  oceans. 

The  level  of  this  component  falls  rapidly  with  increasing  period 
from  the  peak  at  about  16  seconds.  This  portion  of  the  noise  is 
spacially  coherent  and  usually  comes  from  cue  or  two  major 
directions  at  any  given  time.  The  power  level  will  be  similar 
on  horizontal  and  vertical  instruments.  The  power  level  will 
not  be  decreased  by  seismometer  burial  (see  Capon,  1969). 

The  remainder  of  the  noi:  e is  locally  generated  by  atmospheric 
loading.  Turbulent  loading  is  the  most  common  method  of  generation, 
although  occasionally  an  acoustic  gravity  wave  will  significantly  add 
to  the  seismic  noise.  Seismometer  burial  decreases  the  noise 
level  of  the  turbulent  component,  but  not  of  the  gravity  wave 
noise.  For  shallow  seismometers  the  horizontals  will  have  a higher 
noise  level  than  the  verticals. 
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2.  At  the  longer  periods  (30  to  60  seconds),  the 
ocean  generated  microseisms  are  not  significant,  (Savino,  et  al.  , 

1972b;  Capon,  1969;  Sorrells  and  Douse,  197'').  The  major  noi se  source 
is  from  atmospheric  turbulence  static  generation  in  the  vicinity  (apnroxi- 
mately  500  m)  of  the  seismometer.  This  noise  has  the  following  properties 

1.  On  the  vertical  seismometer,  the  minimum 
power  level  is  at  30  to  40-second  period; 

2.  The  power  level  (at  30  seconds)  is  much 
higher  (15  dB)  on  a surface  seismometer 
compared  to  a deep  one  during  turbulent 
atmospheric  conditions,  and  abc  ut  3 dB 
higher  during  calm  conditions.  Power 
levels  vary  by  about  15  dB  with  time; 

3.  Coherence  between  verticals  is  discussed 
by  (Savino,  et  al,  1972b).  The  coherence 

is  on  the  order  of  0.  59  for  a 1 km  separation; 

4.  (Savino,  et  al.  , 1972b,  p.  189).  There  is 
no  significant  correlation  between  the 
horizontal  component  with  the  verticals 
over  the  period  range  of  1 0 to  500  seconds 
due  to  the  noise  being  generated  by  many 
uncorrelated  sources.  (I  believe  that 

in  the  10  to  30  second  band  there  will 
at  t:mes  be  some  coherence  between 
horizontals  and  the  vertical.  ); 

5.  Noise  power  levels  for  OGD  are  similar 
on  horizontal  and  verticals  (Savino,  et  al.  , 

1972b) 
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B.  Earthquake  Interfex  e.ice 

1.  Major  Arrivals 

The  initial , normally  discussed,  arrivals  from 
earthquakes  consist  of  the  following.  No  effort  is  made  to  discuss 
these  here  since  their  properties  are  well  known. 

a*  P and  S body  waves 

b.  Love  waves  (direct  path) 

c.  Rayleigh  waves  (direct  path) 

d.  major  multibounce  body  waves,  PP,  SS,  PPP,  SSS,  etc. 

e*  PL  coupled  leaky  waves. 

These  arrivals  are  all  arrivals  occurring  before  the  Rayleigh 
wave  arrives. 

2.  Multipath  Surface  Waves 
a-  At  20  seconds 

Many  multipath  arrivals  are  observed.  Some  multipath 
arrivals  have  power  levels  comparable  to  the  main  Rayleigh  wave  train. 
Azimuths  of  multipath  arrivals  and  even  the  main  Rayleigh  wave  may 
differ  from  the  great  circle  path  by  as  much  as  40°.  The  multipaths 
arise  because  of  differences  in  phase  velocity  between  different  areas 

of  the  earth  (especially  between  continents  and  oceans).  These  contrasts 
cause  refraction  of  waves. 

In  addition,  if  a wave  hits  a boundary  at  a large  angle  to 
the  normal,  a significant  amount  of  energy  may  be  reflected.  Thus, 
multipaths  can  arise  even  if  source  and  receiver  are  on  the  same 
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continent.  The  best  references  for  these  multipaths  are  Capon(1970) 
and  Capon  and  Everndcn  (1971).  The  interference  of  multipath 

arrivals  is  also  noted  by  Newton  (1973)  and  Boore  and  Toksoz  (1969) 
(See  Table  I). 

b.  At  40  seconds 

The  properties  of  these  multipath  Rayleigh  waves  are 
very  important  for  the  present  work.  The  most  complete  study  of 
their  properties  is  Capon  and  Evernden  (1971).  This  work  is  based 
on  pre-low-pass  filtered  data.  Table  2 gives  a compilation  of 
their  results.  Note  that  the  6dB  bandpass  of  the  filter  used  after  the 
prefilter  was  (.019  - .031  Hz)  or  (32  to  52  second  period). 

From  this  table  it  is  observed  that  the  highest  power 
level  is  always  in  the  first  200  second  window,  or  in  both  the  first 
and  second  windows.  For  the  first  window,  the  azimuth  deviation 
(arrival  angle  difference  from  the  great  circle)  is  0°  with  the  exception 
of  a 10°  deviation  for  Event  ft  17.  Thus,  an  assumption,  that  the 

arrival  direction  of  the  main  40-second  period  wave  is  the  great  circle 
path,  is  justified. 

The  power  levels  of  the  later  arriving  multipaths  is 
generally  down  lOdB  and  usually  much  more  from  the  main  arrival. 

The  multipath  arrivals  do  have  angular  deviations  of  an  order  of  20° 
in  some  cases. 


5 


Table  1:  Measured  Variations  in  Azimuth  of  Arrival  and 

Power  Levels  at  Various  Times  for  Events 
From  Capon  (1970) 
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Table  1:  Measured  Variations  in  Azimuth  of  Arrival  and 

(contd)  Power  Levels  at  Various  Times  for  Events. 
From  Capon  (1970) 
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Table  2 


Measured  Variations  in  Azimuth  of  Arrival  and 
Power  Levels  at  Various  Times  of  the  40-sec 
Period  Rayleigh  Wave  Group  for  10  Events 
From  Capon  and  Evernden  (1971) 
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i L)t(  l 

Power 

l ^ 200 

0 

0 

• 200  — irr) 

1G 

-13 

•100 -Goo 

-5 

-20 

GOO  sort 

-5 

-16 

soo-woo 

— 5 

-10 

1000-1200 

—5 

-22 

1200-HOO 

-i 

-20 

HOO-JCOO 

-§ 

-10 

2 0-200 

0 

0 

200— 1(X3 

0 

—5 

400-G0Q 

0 

-1G 

000  S00 

— I! 

-25 

* 600-1000 

-!! 

—21 

1000-1200 

-II 

-23 

1200  H(0 

HI 

-17 

H00-1G00 

-II 

-21 

3 0-200 

0 

0 “1 

200-100 

0 

-1  J 

400-000 

0 

-21 

CQO-500 

-{ 

-32 

800-1000 

— 30 

-20 

1000-1200 

-30 

-27 

1200-H00 

-§ 

-33 

llOO-lCUO 

-5 

-31 

6 0 200 

0 

0 

200-400 

0 

-9 

400-000 

20 

-9 

G00-800 

20 

-g 

600-1000 

20 

-S 

1000-1200 

-|j 

-11 

1200-1100 

-II 

-13 

HOO-IGOO 

-1! 

-19 

8 0-200 

0 

0 

200-100 

0 

-11 

4OO-GO0 

— II 

-30 

G00-600 

-i 

-31 

800  1000 

-ii 

-30 

1000-1200 

-I! 

-30 

1200-1100 

-II 

-38 

H00-1G00 

-II 

-40 

1G  . 0-200 

0 

0 ""I 

200-400 

12 

0 J 

400  C00 

-II 

vOO-SOO 

25 

800-1000 

-II 

-7 

1000-1200 

-$ 

-13 

1200-1100 

-{ 

-11 

H00-1G00 

—I1 

-8 

similar  power 
in 

two  windows 


8 


Table  2:  Measured  Variations  ir  Azimuth  of  Arrival  and 

(contd.  ) Power  Levels  at  Various  Times  of  the  40-sec 
Period  Rayleigh  Wave  Group  for  10  Events. 
From  C?pon  and  Evernden  (1971) 


Event  No. 
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Dev  utu<nt  • 
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• Time  is 
wave. 

measured  relative  to  onset  of  Rayleigh 

f Azimuthal  deviations 

are  differenecs  between 

measured  azimuth  of  arriva 

1 and  true  azim 

nth  of  cpi- 

center  of  event . 

} Tower  levels  in  decibels  relative  to  power  of 
main  group. 

S l’ower  level  is  below  that  of  background  noise 
find  no  direction  n f arrival  can  be  assigned. 

|'  Power  level  is  above  that  of  background  noise  but 
no  direction  of  arrival  can  be  assigned. 
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The  description  of  Love  wave  multipaths  at  40  second 
period  is  given  by  C?  pon  (1971,  L).  The  Love  wave  multipaths  have 
similar  properties  (azimuth  deviations  and  relative  power)  as  the 

i 

Rayleigh  wave  multipsths. 

3.  Late  Arriving  Earthquake.  Interference 


By  late  arriving  interference,  we  mean  energy 
which  comes  after  the  Rayleigh  wave  and ''simple*  multipaths.  This  inter- 
ference includes 

a)  Surface  waver  going  the  long  way  around 
the  earth  or  going  around  more  than  once, 

b)  Multibounce  body  waves, 

c)  PKPPKP  (also  called  P'P')  waves.  We 
will  be  discussing  the  roughly  40  second 
and  above  period  range. 

Capon  and  Evernden  (1971)  - note  that  the  coda  level  decays 
rapidly  to  more  than  30  dB  below  the  main  Rayleigh  arrival  in 
6-7  minutes  for  some  events.  For  other  events,  the  coda  decreases 
much  more  slowly,  with  coda  only  7 to  13  dB  below  the  main 
Rayleigh  wave  power  for  period  of  5,  000  seconds  or  more.  Rates  of 
decay  are  shown  in  Figure  0,  for  the  first  Z6  minutes. 

By  f-k  determinations  of  azimuth  and  phase  velocity  across  LASA, 
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X AVERAGE  rO'.VcP 
I OH  DATA  ABOVE 
NOISE  LLVEL 
f INDICATES  DATA 
BELOW  NOISE 
LEVEL 


0-33  67-0  113-16  7 20-23! 

33  67  10-13 » 16  7 10  733-207 


TIME  AFTER  ONSET  OF  RAYLEIGH  WAVE  Iran) 


0 3 3 6 7 10  0 13  3 16  7 20  0 23  3 20  7 

DIFFERENCE  IN  ARRIVAL  TIMES  OF  RAYLEIGH  WAVES  (rr  n) 

(b) 


Fig.  0:  Rate  of  Decay  of  Coda  Power  Level  and 

Performance  Characteristic  of  Proposed 
Detection  Method.  From  Capon  and  Evernden  (1V71) 
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Capon  and  Evernden  (1971)  identified  waves  from  the  paths  a),  b),  and 
c)  described  above. 

Savino  et  al.(1972  b)  p 168 

They  give  the  duration  of  coda  from  large  events  (average 
body  wave  magnitude  of  6)  as  approximately  9 hours  at  Ogdensburg,  (OGB), 
but  note  that  the  predominant  periods  are  60-180  seconds.  Thus, 
much  of  this  coda  is  out  of  our  band.  This  type  of  coda  dominates 
the  OGB  record  approximately  10%  of  the  time.  I speculate  that 
with  a high  pass  cutoff  at  60  second  period,  the  fraction  of  time 
that  the  earthquake  coda  predominates  the  atmospheric  noise  would 
be  greatly  reduced. 

Savino  et  ?.l.  (1974  b,  p.  167) 

Following  Alsop  and  Brune  (1965)  it  is  suggested  that  multiply 
reflected  P and  S body  waves  make  up  a large  portion  of  the  late  coda. 
This  agrees  with  Capon  and  Evernden  (1971).  We  note  that  the 
polarization  (ellipticity ) of  body  waves  will  be  much  different  from 
Rayleigh  waves.  If  the  wave  arrives  at  the  seismometer  as  a P wave-> 
it  will  be  predominantly  on  the  vertical;  and  if,  as  an  S wave,  predom- 
inantly on  the  horizontals.  (I  have  not  checked  this  last  statement 
carefully  for  40  second  and  longer  waves,  but  I do  know  of  some 
literature  to  look  at  if  you  desire)  (Fig.  1 gives  OGD  seismograms 
from  seveval  events)* 
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Fig.  1:  Comparison  of  High-Gain  Seismograms  (A),  (B),  and  (C)  at 

OGD  with  a relatively  low-gain  Seismogram  (D)  from  an  instru- 
ment at  Palisades,  New  York,  for  ..he  duration  of  Intermediate 
and  Deep-Focus  Earthquakes.  Epicentral  Data  for  the  Earth- 
quakes in  (A),  (B),  and  (C)  are  given  in  Table  3.  The  Epicentral 
Data  Reported  by  the  USCGS  for  the  Event  in  (D)  are 
Aug.  15,  1963,  17h2 5m06s,  13.  8S,  69.  3W,  Peru,  Bolivia, 
h = 543  km.  The  magnitude  was  estimated  at  7 1/4  by  the 
Seismological  Laboratory,  Pasadena. 

From  Savino  et  al.  1972 


The  general  duration  of  coda  at  40-60  second  period  and  its 
properties  (1.  e.  , types  of  waves)  does  not  seem  to  have  been 
published.  Neither  Savino  et  al.  or  Capon  and  Evernden  have 
addressed  the  very  late  coda  in  this  band  systematically. 

Savino  et  al.  (1972  a) 

They  give  results  from  8 of  the  high  gain  stations.  Detection 
ability  is  discussed,  as  is  noise  levels  at  the  various  stations.  Band 
pass  filtering  is  discussed  as  is  polarization  filtering. 
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2 . 2 Elliptic ity  Values  and  Variations 


Boore  and  Toksoz 

Phase  difference  between  vertical  and  horizontal  components 
were  usually  90°  i 10°  ( p.  339).  This  paper  also  shows  theoret- 
ical phase  between  horizontal  and  verticals  for  multipath  interference. 
Interference  peaks  are  usually  between  20  and  30  seconds. 

Figure  2 shows  ellipticity  values  and  gives  a means  of 
estimating  scatter  as  a function  of  period. 

Newton 

Figure  3 shows  scatter  of  azimuth  of  arrival  data.  f-k  analysis 
with  LASA  gives  about  + 5%  ('90%'  confident"  limited  to  great  circle 
azimuth.  (See  also  Newton,  p.  c3).  Using  amplitude  ratios  of  horiz- 
ontal power  of  f-k  gives  errors  in  azimuth  of  about  8°  both  positive 
and  negative  This  work  was  done,  however,  on  only  a few  carefully 
selected  earthquakes,  and  applied  to  the  main  part  of  the  Rayleigh 
wave  train. 

Higher  mode  Rayleigh  waves  were  also  studied.  These  waves 
have  higher  group  velocity  than  the  fundamental  mode  ar.d  arrive 
somewhat  before  the  fundamental  mode.  The  eFipticity  for  these  higher 
modes  differs  from  the  ellipticity  of  the  fundamental  mode.  (Figure  4, 
a,  b,  and  c) 

On  Figure  5,  I have  compared  the  Boore-Toksoz  results  with 
Newton.  Note  that  they  differ. 
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ELUPTXITY  »M/V) 


Fig.  2:  Measured  Ellipticity  for  Different  Events  and 

Different  Sites.  The  Solid  Line  is  the 
Theoretical  Ellipticity  Corresponding  to 
Model  Ml.  From  Boore  and  Toksoz. 
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Fig.  3:  Results  of  the  f-k  Analysis  of  Signals  from  the  Tonga  Isla 
and  Greenland  Sea  Earthquakes.  (a)  Normalized  Amplitud 

Spectra.  (b)  Phase  Velocities.  (c)  Propagation  Direction 
(From  Newton,  1973,  p.  109) 
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Fig.  4a:  Composites  of  Ellipticities  Observed  at  the  .CASA  Subarrays 

in  Group  I.  From  Newton  (1973) 
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in  Group  II.  From  Newton (1973) 
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Fig.  4c:  Higher  Mode  Ellipticities  from  Observations  at  LASA  Subarrays  of  the 

Tonga  Island  Earthquake  Signals.  Open  Symools  Denote  Reciprocal 
Ellipticity  (i.  e.  , Vertical  over  Radial  Spectral  Amplitudes) 

Note:  These  will  usually  arrive  before  the  Fundamental. 


Approximate  range  of  Ellipticity 
values:  of  Boore-Toksoz  compar< 


o , j Period 

Comparison  of  Boore-Toksoz  Results  with  Newton  Res 


and  the  arrival  time  of  a given  frequency  is  given  by 

t«.  ■ y / U(J) 

where  (J  (wl)  is  the  group  velocity,  given  by 

For  the  40  to  60  second  period  range  (J  ( S)  is  fairly 
constant.  Thus,  energy  arrives  over  a fairly  short  period  of  time. 
Figures  6 and  7 show  LASA  LPZ  seismograms  after  filtering  by 
the  filter  shown  in  Figure  8.  For  the  event  of  Figure  6,  the  energy 
is  contained  in  a 200  second  window.  Thus,  a 200  second  window,  as 
presently  being  used  by  OAS,  is  sufficient  for  this  event.  However, 
it  remains  to  develop  a method  to  select  tha  proper  time  interval, 
Overlapping  windows,  starting  every  30  seconds  would 
probably  assure  one  window  covered  the  arrival.  Perhaps,  the 
criteria  of  the  window  with  the  most  energy  in  the  explosive  event 
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Event 


The  Prefiltered  Wave  Forms  for  Event  1 Added 
Artifically  to  Those  of  Event  3. 

From  Capon  and  Evernden  (1971) 


06  41  40 


08  10  40 


06  11  40 


1 I 

09:00  40 


(b) 


Fig.  7:  Prefiltered  Wave  Forms  of  Event  19  of  October 
15,  1967,  near  the  Coast  of  Nicaragua,  Snowing 
Initial  Rayleigh-Wave  Group  Arrival  and  Multiply 
Reflected  P Waves.  From  Capon  and  Evernder.  (1971) 


TIMC  (mm) 


Fig.  8:  Frequency  Response  and  Impulse  Response  of 
of  Fifth  Order  Low-Pass  Chebyshev  Filter. 
From  Capon  and  Evernden  (1971)  Pro  Low- 
Pass  Filter. 
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azimuth  direction  (cn  the  rotated  horizontal)  could  be  used  to 
determine,  which  of  the  windows  actually  contains  the  arrival. 

If  the  window  used  does  not  contain  most  of  the  energy  in 
the  arrival,  I would  expect  an  error  in  the  measured  phase  spectrum. 
Figure  7 shows  an  event  (#19)  with  a Rayleigh  wave  of  300  seconds 
length.  From  Table  2,  we  see  that,  for  three  (Events  3,  16,  and  19) 
of  the  ten  studied,  the  power  levels  were  similar  for  the  first  two 
200  second  windows.  For  these  events,  a longer  window  should  be 
used. 

For  the  20-40  second  band,  the  group  velocity  increases 
with  period.  This  portion  of  the  Rayleigh  wave  train  is  highly 
dispersed  and  typically  arrives  over  a 400-600  second  period  (the 
arrivals  are  extended  even  more  in  time  by  multipaths).  Thus,  lor 
this  period  range  no  single  200  second  window  will  contain  all 
frequencies.  Results  of  Capon  and  Evernden  (1971,  p.  819)  and 
Newton  (1973)  demonstrate  this  point.  Also  see  Table  I and  records 
shown  in  Figures  9-14. 

I have  not  presently  developed  this  approach,  but  rather  than 
using  spectral  analysis,  the  bank  of  bandpass  filters  approach  with 
cross  correlations  between  components  to  separate  the  Rayleigh  waves 
of  two  events  might  be  considered. 

For  the  LASA  (25  second  Peak!  seismometers,  20  second 
period  energy  leaked  through  the  digital  filters  into  the  40  second 
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Fig.  1 1: 


The  Long-Period  Wave  Forms  for  Event  19, 
the  October  15,  1967,  Event  near  the  Coast 
of  Nicaragua.  From  Capon  and  Evernden  (1971) 
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Fig.  12:  Narrow- bandpass  Filter  Traces  cf  the  Vertical 

Component,  Mid-Atlantic  Ridge  .Earthquake. 
From  Newton,  1973. 
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Fig.  13:  Narrow-bandpass  Filter  Traces  of  the  Radial  Component 
Seismogram  Containing  the  Fundamenta1  Mode  Rayleigh 
Wave  from  the  Mid-Atlantic  Ridge  Earthquake. 
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Periods  (sec)  at  peak  of  Gaussian  filter  transfer 
functions . 
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Fig.  14:  Narrow-bandpass  Filter  Traces  of  the  Transverse  Component 
Seismogram  Containing  the  Fundamental  Mode  Rayleigh  Wave 
from  the  Mid-Atlantic  Ridge  Earthquake. 
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spectral  estimates  used  for  the  high  rssolution  f-k  analysis.  This 
occurs  because  the  20  second  power  level  is  much  above  the  40 
second  power  on  the  LASA  seismometer  outputs.  To  eliminate  this 
problem,  Capon  and  Evernden  (1971)  used  a low  pass  Chebyshev 
filter.  (See  Figure  8).  The  spectral  leakage  problem  will  probably 
not  be  a problem  with  our  high  gain,  longer  period  seismometers. 
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3.  2 On  the  Ellipticity  Filter 


Newton  (1973)  (copy  of  appendix  enclosed)  tried  a method  of 

determining  the  arrival  azimuth  based  on  the  assumption  of  a 90°' 

phase  shift  between  the  radial  and  vertical  signals.  The  method  was  to 

find  the  direction  at  which  the  phase  shift  was  9C°.  The  method  worked 

for  simulated  data.  However,  it  did  not  do  well  on  real  signals. 

o 

Newton  attributed  the  failure  to  an  extreme  sensitivity  to  the  90 

phase  shift  assumption.  He  points  out  that,  if  the  ellipticity  is 

changing  with  frequency,  the  "ellipticity"  of  a finit  bandwidth  signal 

o 

will  not  have  a 90  phase  shift.  He  also  did  an  experiment  in  which 
artificial  data  was  used.  When  the  radial  data  was  delayed  by 
one  second,  the  phase  shifts  which  resulted  led  to  large  errors  in 
the  azimuth  determined  by  the  algorithm. 

The  difficulty  Newton  had  in  getting  the  azimuth  for  a single 
Rayleigh  wave  may  indicate  that  we  will  have  difficulty  in  the  case  of 
the  arrival  of  two  Rayleigh  waves.  If  we  do  have  difficulties,  one 
thing  that  may  help  is  to  assume  that  the  explosion  signal  is  coming 
along  the  great  circle  direction.  (This  assumption  is  valid  for  most  events 
at  40  second  periods  and  longer;  this  was  discussed  in  a previous 
section  of  this  report)  The  five  other  quantities  (complex  amplitudes 
of  the  two  signals  and  azimuth  of  interfering  events)  could  be  found  by 
least  squares.  This  procecu  *e  would  perhaps  be  more  robust  than 
the  solution  for  all  6 quantities.  This  approach  is  supported  by 
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Newton  (p.  93);  he  states  that  the  most  correct  estimate  of  the  radial 
power  was  obtained  by  using  the  great  circle  azimuth. 

It  appears  that  multiple  reflected  P and  S body  waves  will 
also  cause  difficulty.  It  may  be  worth  while  to  consider  the  effect 
of  these  waves  on  the  filter. 
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4.  0 General  Comment 


I believe  that  there  are  limits  to  the  signal  - to  interference 
gain  - we  can  get  with  a single  seismometer.  In  the  course  of  our 
work  on  the  single  seismometer  processing,  I believe  we  should 
assemble  data  on  signal  properties  and  the  three-component  coherencies 
of  signals.  If  we  can,  we  should  consider  processing  mrfhods  for 
small  arrays,  with  perhaps  4 high  gain , 3 -component  seismometers. 

If  the  40  second  band  pass  is  useful,  then  I am  sure  arrays  will 
eventually  be  used. 

I think  that  40  second  periods  may  be  of  great  use  in  source 
estimation.  This  is  because  lateral  refractions  are  not  important 
at  this  frequency.  Thus,  geometrical  spreading  and  focusing  (McGarr, 
1969;  McGarr,  1972;  Aki  et  al.  1972)  which  is  so  severe  at  20  second 
period  will  not  be  as  important  at  40  seconds.  Therefore,  amplitude 
estimates  at  40  seconds  should  be  more  reliable. 

I have  seme  information  at  hand  on  explosion  spectra  but 
have  not  included  it  at  this  time. 
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Appendix 

(From  Newton  1973) 
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Appendix  E 

NULL  METHOD  FOR  DETERMINING  RAYLEIGH  WAVE 
PROPAGATION  DIRECTIONS 

The  fundamental  and  higher  modes  of  Rayleigh  waves  are 
characterized  by  the  elliptical  particle  motion  where  the 
vertical  and  horizontal  components  have  a constant  phase  dif- 
ference of  ir/2.  If  the  motion  is  retrograde  the  radial 
phase  leads  and  for  prograde  motion  the  vertical  phase  leads, 
and  the  transition  from  one  to  the  other  occurs  at  a node  of 
either  of  the  components. 

If  w (t)  denotes  a vertical  component  seismogram  then 
the  analytical  signal  is  the  complex  representation  of  w(t), 
defined  by  S(t)  = w(t)-iw(t),  where  w(t)  is  the  Hilbert 
transform  of  w(t).  w(t)  is  referred  to  as  the  "quadrature" 
signal  and  w(t)  as  the  "cophase"  signal,  (see  Bracewell, 

1965) . The  cophase  and  quadrature  signals  are  also  charac- 
terized by  a ir/2  phase  difference,  and  the  quadrature  phase 
leads . 

Now  consider  the  horizontal  component  seismograms  h! 
and  h 2,  which  are  recorded  by  the  orthogonally  oriented 
instruments.  If  6'  is  the  angle  between  the  installed  azimuth 
for  hi  and  the  propagation  direction,  0 , for  a Rayleigh 
wave,  then  the  radial,  u(t),  and  transverse,  v(t),  seismo- 
grams are 


r 

“■  *■ 

u 

COS0' 

sinO  1 

hi 

= 1 

• 

V 

1 -sin0 ' 

cos  0 1 
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Assuming  that  only  Rayleigh  motion  is  present  on  u and  w 
then 


u-iu  U 

7T  = w 


i(4>  -<J>  ) 
Yu  Yw 


I U . ±iir/2 

lwle 


*ilSH 


Therefore 


w-iw 


Reai\Hii2  = = 0 


v 2+w  2 


which  gives  the  null  relation 


uw+uw  = 0. 


Making  the  substitutions 


= hj  cos  0 '+h  2 sin  0 ' 


h4  cos  0 '+h  z sin  0 ' , 


we  can  solve  for  01: 


0 1 = -tan" 


(h  w+h  w) 

i i 
(h  2W+h  2W) 


This  algorithm  was  tested  using  synthetic  signals,  the 
radial  component  having  the  same  spectrum  and  group  delays 
as  the  vertical  but  a constant  phase  difference  of  tt/2,  and 
the  transverse  component  was  sinusoidal  with  slowly  changing 
period  and  with  maximum  amplitude  approximately  the  same  as 
the  other  components.  Except  for  the  leading  and  trailing 
edges  of  the  Rayleigh  wave  window  the  computed  0'  was  zero. 
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When  the  vertical  component  was  shifted  by  one  second  with 
respect  to  the  horizontal  components  O'  was  seldom  zero  and 
only  for  a short  window  where  the  signal  amplitude  was  the 
largest  did  O'  usually  fall  within  10  degrees  of  zero.  This 
means  that  digitized  seismograms  must  have  nearly  perfect 
synchronization  between  componrnts,  i.e.,  instrument  group 
delays  must  be  uniform  and  sampling  must  be  synchronous. 

Another  pitfall  for  this  algorithm  occurs  when  the  el- 
lipticity  is  changing.  This  effect  can  be  seen  graphically 
by  drawing  two  vectors  having  a given  amplitude  ratio  and 
phase  separation  and  then  drawing  two  vectors  with  phases 
ir/2  greater  than  the  first  two  and  with  a different  ampli- 
tude ratio.  Although  the  tt/2  phase  difference  exists 
between  the  pairs,  it  does  not  exist  between  the  resultant 
vectors.  It  was  thought  that  this  problem  should  vanish 
when  band  filtering  is  used  because  the  ellipticities  for 
realistic  earth  structures  are  slowly  varying  with  frequency. 
However,  the  complete  failure  for  the  method  when  applied  to 
high  signal-to-noise  seismograms  prompted  a test  of  this 
assumption. 

Using  the  cophase  and  quadrature  signals  of  a vertical 
component  seismogram  as  the  vertical  and  radial  test  signal 
components,  it  was  found  that  uw  + uw  did  indeed  vanish. 

When  broad-band  Gaussian-shaped  filters  were  applied  to  the 
two  comp'  .ents  using  slightly  different  center  frequencies, 
uw  + uw  did  not  vanish.  It  was  concluded  that  slight  depar- 
tures from  a tr/2  phase  difference  between  components,  such 


40 


as  would  appear  with  even  small  amounts  of  noise  and  when 
the  ellipticity  is  frequency  dependent,  arc  sufficient  to 
render  useless  the  null  method  for  determining  Rayleigh  v 
propagation  directions. 
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